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Abstract 

Brachydactyly type A2 (BDA2, MIM 112600) is characterized by the deviation and shortening of the middle phalange of the 
index finger and the second toe. Using genome-wide linkage analysis in a Chinese BDA2 family, we mapped the maximum 
candidate interval of BDA2 to a ~1.5 IVlb region between D20S194 and D20S115 within chromosome 20pl2.3 and found 
that the pairwise logarithm of the odds score was highest for marker D20S156 (Zmax = 6.09 at 9 = 0). Based on functional 
and positional perspectives, the bone morphogenetic protein 2 (BMP2} gene was identified as the causal gene for BDA2 in 
this region, even though no point mutation was detected in BMP2. Through further investigation, we identified a 4,671 bp 
(Chr20: 6,809,218-6,813,888) genomic duplication downstream of BMP2. This duplication was located within the linked 
region, co-segregated with the BDA2 phenotype in this family, and was not found in the unaffected family members and 
the unrelated control individuals. Compared with the previously reported duplications, the duplication in this family has a 
different breakpoint flanked by the microhomologous sequence CATCA and a slightly different length. Some other 
microhomologous nucleotides were also found in the duplicated region. In summary, our findings support the conclusions 
that BMP2 is the causing gene for BDA2, that the genomic location corresponding to the duplication region is prone to 
structural changes associated with malformation of the digits, and that this tendency is probably caused by the abundance 
of microhomologous sequences in the region. 
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introduction 

The term brachydactyly (BD) is derived from the ancient Greek 
(brachy-: short; dactylos: digit). This term indicates the shortening 
of digits due to abnormal development of tlie phalanges and/or 
the metacarpals and may also be accompanied by other hand 
malformations, such as syndactyly, Polydactyly, reduction defects, 
or symphalangism [1]. One of the most commonly used 
classifications of brachydactyly based on anatomical grounds was 
provided by Bell [2] and further elaborated by Temtamy and 
McKusick [3] . BD can occur either as an isolated malformation or 
as part of a complex malformation syndrome. The isolated BDs 
are classified clinically into BDA-E and three other subgroups [3] . 
BD types have certain overlapping features, including hypoplasia/ 
aplasia of phalanges and abnormal interdigital joint formation, 
suggesting that the formations of phalanges and joints are linked 
on a developmental and molecular basis. 



In recent years, most of the isolated brachydactylies have been 
characterized at the molecular level, and these studies have 
provided new insights into the mechanisms of joint and digit 
development [4]. BDAl (MIM 112500) is characterized by the 
shortening or absence of the middle phalanges and is associated 
with mutations in the Indian hedgehog gene {IHH) or chromo- 
some 5pl3.3-pl3.2 [5,6,7,8,9]. BDBl-2 (MIM 113000; MIM 
61 1377) is characterized by terminal deficiency of the fingers and 
toes and mutations that have been identified in an orphan receptor 
tyrosine kinase ROR2 or in NOGGIN [NOG) [10,1 1,12,13]. BDC 
(MIM 1131 00) can present with a range of anomalies (from 
shortening to hyperphalangy of the middle phalanges of the index, 
middle, and littie fingers) and is caused by mutations affecting 
gTowth/differentiation factor 5 {GDF5) [14]. Two mutations of 
H0XD13 are associated with distinctive limb phenotypes in which 
brachydactylies of specific metacarpals, metatarsals, and phalan- 
geal bones are the most constant features and exhibit overlap with 
brachydactyly types D (MIM 113200) and E (MIM 113300) [15]. 
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BDA2 (MIM 1 12600) was described first by Mohr and Wiredt 
in a large Norwegian kindred of Danish descent [16]. Patients 
afflicted with BDA2 display hypoplasia/aplasia of the middle 
phalanx of the second and, sometimes, fifth fingers. Common to 
all BDA2 families are the deviation and shortening of the middle 
phalange of the index finger and/or the second toe. BDA2 is 
genetically heterogeneous. Lehmann et al. found that BDA2 is 
caused by missense mutations in bone morphogenetic receptor 1 b 
(BMPRIB), a type I transmembrane serine-threonine kinase, and 
that three different missense mutations exist in three German 
families [17,18]. The second BDA2 -causing gene, GDF5, was 
identified in the descendants of the original BDA2 family by 
Seemann et al. and Kjaer et al. [19,20,2 1] . Recendy, we discovered 
a Chinese BDA2 family, but no mutation was found either in the 
BMPRIB or GDF5 loci in the genetic analysis, suggesting 
additional genetic heterogeneity of this trait. Therefore, we 
performed a genome-wide scan using 403 microsatellite markers 
in our BDA2 family. 

More recendy, duplications downstream oi BMP2 (NM 001200) 
were identified in two studies; one study comprised samples from 
two European families (~5.9 kb and ~5.5 kb), and the other study- 
comprised samples from a Chinese family (~4.6 kb) [22,23]. In 
our current study, we found a 4,67 1 bp duplication downstream of 
BMP2 (Chr20: 6,809,218-6,813,888) associated with BDA2 in our 
Chinese family that partially overlaps with the previously reported 
duplications but has a different breakpoint flanked by the 
microhomologous sequence GATCA. Our finding supports the 
hypothesis that the genomic location corresponding to the 
duplication region is prone to structural changes associated with 
malformation of the digits and that this tendency is probably 
caused by the abundance of microhomologous sequences in the 
region. 

Results 

Clinical phenotypes 

In this six-generation family, 19 affected and 10 unaffected 
individuals were clinically examined, and a subset of the affected 
individuals was selected for radiography. AH of the affected 
individuals had the t)'pical phenotypes of BDA2 characteristic of 
medially deviated and shortened index fingers and second toes 
with abnormal interdigital joint formation. Most of the affected 
adult individuals exhibited short stature. Triangular-shaped 
middle phalanges were observed in all of the affected individuals, 
as demonstrated in the photograph and X-rays of two affected 
adult (Figure 1). No sex-correlated characteristics were observed in 
the family. 

The common characteristic of BDA2 caused by variants in 
different genes is that almost all of the patients showed shortened 
index fingers. The phenotypic differences consisted of some 
abnormal symptoms in other phalanges or metacarpals in some 
patients and a few severely affected individuals, indicating overlap 
with another type of brachydactyly, symphalangism, or syndactyly. 
In addition, few mutation carriers appeared to be clinically normal 
(Table SI). No variant-specific phenotypes were observed in the 
affected individuals. 

Genetic analyses of BMPRIB and GDF5 

The sequence analyses of the known BDA2 causing genes, 
BMPRIB and GDF5, in four affected famUy members (111:2, IV:4, 
IV: 12, and V:12 in the family) revealed no mutations. In the two- 
point linkage analysis, the LOD scores were calculated for 
microsatellite markers in the BMPRIB and GDF5 loci, but none 
of the LOD scores exceeded 2.0 (at 9 = 0) (data not shown). 



Genome-wide scan and linkage analyses 

A two-point LOD score for 21 markers on chromosome 
20pl2.3 for the pedigree are summarized in Table 1. Seventeen 
adjacent markers provided a positive LOD score (9 = 0), and the 
locus for BDA2 was mapped between markers D20S867 and 
D20S851 witiiin a region of -5.2 Mb, based on the UCSC 
Human Genome (Mar 2006). The pairwise logarithm of the odds 
scores was highest for marker D20S156 (Zmax = 6.09 at 9 = 0). 

The subsequent haplot^pes were constructed using the 
CYRILLIC 2.1 software to define the interval of the linked region 
(Figure 2). The markers (19 total) used are hsted in Table 1, and 
the haplotypes were checked through visual inspection. The 
recombination events between the BDA2 phenotype and the 
markers spanning the region of interest defined the smallest 
cosegregating region that included critical meiotic recombinants in 
the pedigree. Careful examination of the haplotypes confirmed 
that the disease-associated alleles co-segregated with the pheno- 
type of BDA2 in the pedigree. A recombination event in individual 
11:2 placed the disease locus distal at D20S194 because the affected 
individual did not inherit the disease-linked D20S194 alleles 
observed in the family. The other recombination event in 
individual IV: 13 placed the disease locus close to D20S115 
because this unaffected male shared the disease-linked alleles of 
D20S1 15 with his affected mother. Hence, the maximal interval of 
linkage with the BDA2 phenotype is bordered by D20S194 
(telomeric) and D20S 1 1 5 (centromeric) within a region of 
~1.5 Mb. Inspection of this interval revealed only one gene, 
BMP2. 

Mutation and duplication analyses of BMP2 

The sequence analysis of BMP2 (i.e., all exons, introns, 
promoter region, 5'UTR, 3'UTR, and evolutionary conserved 
regions) showed no mutation. 

Using the Chromosome 20 array, a ~5.0 kb duplication 
(Chr20: 6,809,000-6,814,000) was identified in an affected famUy 
member (IV: 12) downstream oi BMP2 (-110 kb; Figures 3 A and 
4A). 

Utilizing Q;PCR with three primer pairs (P2-P4) within the 
duplicated region and two primers (P 1 and P5) within the flanking 
sequences (Figures 3B and 4B), this duplication was confirmed in 
eight affected famUy members (IV:3, rV:4, IV:7, IV:12, IV:15, 
V:9, V:10, and VI: 1) and was not found in eight unaffected family 
members (111:7, IV:1, IV:6, IV:8, IV:13, V:3, V:7, and V:8). 

Using primers F7 and R7, a 655 bp PCR fragment was only 
produced in the patients (Figures 4B and 5). Moreover, through 
Q;PCR and breakpoint-PCR, this duplication was not found in 
150 healthy control individuals excluding a common copy number 
\ ariant (CNV). Through the direct sequencing of this 655 bp 
fragment in idght affected family members (IV: 3: File SI; rV:4: 
File S2; IV:7: File S3; IV:12: File S4; IV:15: File S5; V:9: File S6; 
V:10: FUe S7; and VI:1: File S8), we identified a 4,671 bp 
duplication (Chr20: 6,809,218-6,813,888) downstream oi BMP2 
and found that this duplication is flanked by the microhomologous 
sequences GATCA on die 5' end (Chr20: 6,809,213-6,809,217) 
and 3' end (Chr20: 6,813,884-6,813,888). Furtiiermore, in 
addition to the rs6085706 SNP, another microhomology (52 bp) 
could also be identified at the 5' (Chr20: 6,809,160-6,809,211) 
and 3' (Chr20: 6,813,831-6,813,882) flanking region of the 
4,671 bp duplicated sequence. In addition, the microhomologous 
nucleotides GTGACC (7 bp) [22] and TTGCAGTGAGC (1 1 bp) 
[23] reported previously were also observed downstream of 
GATCA. The comparison with the previously reported duplica- 
tions showed that the duplication identified in this study does have 
a different breakpoint, overlaps with the other three duplications. 
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Figure 1. Phenotype of BDA2. A and B: Short and medially deviated index finger (above) and second toe (below). C and D: Radiograph of the 
triangular-shaped middle phalanx of the index finger (above) and the second toe (below). 
doi:1 0.1 371 /journal.pone.0094201 .gOOl 



and difTers in its 5' and 3' ends (Figure 6). Based on the above- 
described findings, we conclude tliat the fragment was duplicated 
and is associated with BDA2. 

Discussion 

The family described in the current study has typical 
characteristics of BDA2 [1,24,25]. We excluded BMPRIB and 
GDF5, which were found to cause BDA2 in the previous studies 
[17,18,19,20,21], and mapped BDA2 to the locus within 
chromosome 20pl2.3 because a strong linkage was observed 
between the BDA2 locus and the microsatellite markers in this 
chromosome region. The recombination events in individuals 11:2 
and IV: 13 of the pedigree reduced the interval to a region of 
-1.5 Mb between D20S194 and D20S115. 

Two factors enabled us to identify BMP2 as the BDA2-causing 
gene in chromosome 20pl2.3. First, the results of genome-wide 
scan showed that BMP2 is the only known gene in the candidate 
interval, i.e., this gene is located at 6.697-6.709 Mb within the 
candidate interval (6.65-7.61 Mb on chromosome 20pl2.3). It has 
been estimated that ~25% of the human genome consists of long 
regions that contain no protein-coding sequences [26]. Some of 
these long regions contain regulatory sequences that act as large 
distances to control the expression of neighboring genes, and the 
sequence elements play critically important and conserved 
biological roles [27,28,29,30]. Disruption of the long-range control 
of gene expression may cause disease [31]. BMP2 was embedded 
in a long region with no known genes ~0.6 Mb upstream and 
~ 1 . 1 Mb downstream that was linked with the phenotype of 
BDA2 in the pedigree. It has been reported that the transcription 
oiBMP2 is regulated by an enhancer region located 156.3 kb from 
the promoter of BMP2 and that its noncoding sequence variants 
may influence BMP2 transcription and thus may contribute to 
bone mass and osteoporosis [32,33]. 

Second, BA-IP2 has a functional role in the BMP signaling 
pathway similar to that of the previously identified genes, i.e., 
BMPRIB and GDF5, that cause BDA2 disease 
[17,18,19,20,21,34]. The proteins of the transforming growth 



factor fi (TGF-fi) superfamUy are involved in the process of bone 
and joint development, cell proliferation, and differentiation. The 
superfamUy encompasses extracellular ligands, including bone 
morphogenetic proteins (BMPs) and growth and differentiation 
factors (GDFs) [35]. The BMP type I receptors (BMPRIA and 
BMPRIB) and one type II receptor (BMPR2) are the receptors of 
BMPs and are needed to form a functional complex to initiate 
further signaling events [36,37]. Specific mutations in GDF5 that 
result in a selective inactivation of the binding to BMPRIB are a 
cause of BDA2 [19,20], similarly to dominant negative mutations 
in the GDF5 high-aflfmity receptor BMPRIB itself [17,18]. BDA2 
has also been shown to be caused by mutations in GDFS that 
interfere with cleavage of the GDF5 prodomain from the mature 
peptide [21]. BMP2 activates the BMP pathway by binding to the 
receptor complex [38] and plays an essential role in cartilage 
differentiation and bone formation, similarly to GDF5 and 
BMPRIB [17,19,39,40,41]. It is known that the global loss of 
BMP2 in an animal results in an embryonic lethality and that 
BMP2 plays a center role in skeletal development, such as 
chondrogenesis, cartilage and bone development, the formation of 
bone mineral density, and limb patterning [39,42,43]. Therefore, 
our linkage results and the reported functional studies strongly 
support the conclusion that BMP2 is the causal gene of the BDA2 
phenotype presented in the family analyzed in the present study. 

Through further investigation, we identified a 4,67 1 bp (Chr20: 
6,809,218-6,813,888) genomic duplication downstream oi BMP2. 
This duplication is located within the linked region and was 
cosegregated with the BDA2 phenotype in this family but not the 
unaffected family members nor the unrelated control individuals. 
Recently, Dathe et al. reported two microduplications (~5.5 and 
-5.9 kb) (Chr20: 6,808,477-6,814,024 and Chr20: 6,808,129- 
6,814,024) in a noncoding sequence ~110kb downstream of 
BMP2 in two European BDA2 families. Their study found that the 
duplication contains highly conserved sequences that constitute a 
cis-regulatory element regulating BMP2 expression in the limb and 
that this regulating element duplication can be considered a 
mutational mechanism for BDA2 [22]. More recently, Su et al. 
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Figure 2. Pedigree structure and hapiotypes of the family. The marker order was determined from the Marshfield map and the UCSC Human 
Genome database (Mar 2006). The open symbols indicate the unaffected individuals, the blackened symbols indicate the affected individuals, the 
squares indicate the males, and the circles indicate the females. The blackened bars indicate the chromosome region shared by the affected members 
of the pedigree. 
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identified a duplication of -4.6 kb (Chr2(): 6,809,382-6,814,044) 
in the same region in a Chinese family affected with BDA2 [23]. 
Our work shows a 4,671 bp duplicated sequence in the affected 
individuals and a 4,507 bp overlapping region (Chr20: 6,809,382- 
6,813,888) with the previous studies. A microhomology GATCA is 
located at the 5' and 3' flanking region of this duplicated 
sequence. The generation of the nonrecurring rearrangement 
agrees well with the replication-based model called fork stalling 
and template switching proposed by Lee et al. [44] . According to 
this hypothesis, the microhomology identified in the present study 
may confuse the DNA replication machinery, causing the 
replication fork to stall or pause, and lead to a switch of the 
lagging strand to another active replication fork in physical 
proximity before resuming replication of the original DNA 
template, thereby resulting in duplication formation. Furthermore, 
a longer microhomology (52 bp) was identified upstream of 
GATCA, and microhomologies GTGACC (7 bp) [22] and 
TTGCAGTGAGC (11 bp) [23], which are associated with the 
duplications identified in two other BDA2 families, were also 
observed downstream of GATCA in our study. Therefore, these 
findings support the conclusion that this region of the genome is 
prone to structural changes associated with malformation of the 
digits and that this tendency is probably caused by the abundance 
of microhomologous sequences in the region. 



In a functional study, Dathe et al. observed that a 5.5 kb tandem 
duplication in BMP2 may increase the expression of BMP2 in the 
Umb of transgenic mouse. They speculated that the overexpression 
of BMP2 may result in a deregulation of the fine-tuned BMP 
pathway by increasing the BMPRIA signal and in a relative 
decrease of BMPRIB signaling, which is compatible with the 
previously proposed molecular pathology of BDA2 [22]. Never- 
theless, Su et al. revealed reduced activity of conserved 2.1 kb 
sequences in the duplication region in osteosarcoma U-20S and 
HeLa cells through luciferase activity assays [23]. The discrepant 
expression caused by the duplication downstream of the BMP2 
gene in vivo and in vitro may be due to the complex regulation of 
BMP2 and the other malformations ofjserved in addition to BDA2 
in some of the family members. It has been found that BMP2 
contains two promoters and that the transcription-initiation site is 
not common between the human and mouse BMP2 genes, 
although the sequence of the BMP2 gene is well conserved in the 
promoter region between the two species [45]. In the 2T3 
osteoblast cell line, a proximal promoter fragment can be 
autoregulated by BMP2 [46]. The studies of a ~2.7 kb BMP2 
promoter fragment have suggested that NF-kB and Gli2 can bind 
separately to it and thus drive BMP2 expression in different cells 
[47,48]. Moreover, the 3'-UTR of BMP2 contains essential 
regulatory elements, and many serial regulatory elements are 
located within the large conserved non-coding regions of BMP2 
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the affected individuals and in a ratio of 1 .5 relative to the two copies of the healthy control. The localization of the Q-PCR amplicons is illustrated in 
Figures 4A and 4B. 
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[49,50]. Therefore, the expres.sioii of BMP2 regulated by the 
duplication sequences may depend on different experimental 
conditions and/ or interaction with other cis-regulatory elements. 

BDA2 has genetic heterogeneity. The previous and current 
studies on BDA2 families show that three missense mutations in 
BMPRIB, two missense mutations in GDF5, and three duplication 
elements approximately 1 1 0 kb downstream of BMP2 can 
separately cause the phenotype of BDA2 [17,18,19,20,21,22,23]. 
Despite the complex molecular mechanisms, these findings suggest 
that these three genes are all involved in the BMP signaling 
pathway and that the cause of BDA2 is a shift in the signaling 
intensity from BMPRIB to BMPRIA [2,5]. BMP signahng has 



functions in the growth plate at various stages of differentiation for 
cartilage and bone formation because it can interact with the 
BDAl-causing gene IHH, some TGF-fi proteins, and other 
extracellular and nuclear factors in endochondral ossification 
and thereby form a molecular regulatory network 
[42,51,52,53,54]. Under different conditions, BMP2 can mediate 
BMP signaling through the canonical SMAD pathway to induce 
SMAD target genes, through a mitogen-activated protein kinase 
(MAPK) pathway to result in the induction of alkaline phosphatase 
(ALP) activity, or through a protein kinase C (PKC) signaling 
pathway to increase Bax/Bcl-2 caspase activity [37,38,55]. 
Therefore, there is no simple interpretation for BMP2 function 
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the breal<point region. The numbers above the BMP2 denote the genomic position. Tel: telomere; Cen: centromere. B. Structure of duplicated region 
in patients. The duplicated region is between Chr20: 6,809,218 and 6,813,888. PCR using primers F7 and R7 amplified the 655 bp sequences spanning 
the breakpoint. C. Sequencing results of the breakpoint PCR product. The nucleotides of microhomology A (5 bp) are marked by a green box, and the 
nucleotides of upstream microhomology B (52 bp), including the rs6085706 SNP, are underlined. In addition, the microhomologous nucleotides 
GTGACC (7 bp) [22] and TTGCAGTGAGC (1 1 bp) [23] reported previously were also observed downstream of microhomology A (data not shown). 
doi:1 0.1 371/journal.pone.0094201 .g004 



in the phenotypes of BDA2, and further efforts to clarify the 
regulation and the role of BMP2 in the BMP signaling pathway 
and molecular network are warranted. 

In summary, we identified a 4,671 bp duplication with a novel 
breakpoint flanked by the microhomology GATCA downstream of 
BMP2 in a Chinese family affected with BDA2. Our findings 
support the conclusion that the genomic location that corresponds 
to the duplication region is prone to structural changes associated 
with malformation of the digits and that this tendency is probably 
caused by the abundance of microhomologous sequences in the 
region. Additionally, our study further exemplifies the peculiar 
nature of the balance of expression of the BMP2 gene during 
development, and the genomic location corresponding to the 
duplication region is critical for the function of this gene. 
IMoreover, our results indicate that such duplications may be 
considered a regulatory mutational mechanism. Therefore, similar 
duplications may conceivably be implemented into bioinformatic 
search tools to detect variants as an additional method to associate 
a genetic mutation with a given phenotype. 



Materials and Methods 

Subjects 

The study was approved by the genetic research ethics 
committees of Shanghai Jiaotong University. A Chinese family 
with BDA2 was recruited for the study (see Figure 2). This family 
was from Anhui province in eastern China and does not have any 
distant relatives. The informed consent was written and obtained 
from all participants. 

Genetic analyses of BMPRIB and GDF5 

Samples of peripheral blood DNA were collected from all 
available family members, and the DNA was isolated using 
standard procedures. We amplified the coding region, the spUce 
junctions, and the promoter sequence of the BMPRIB and GDF5 
genes using standard protocol polymerase chain reaction (PCR). 
The PCR products were analyzed using 1.5% agarose gel 
electrophoresis. The purified PCR products were bidirectionally 
sequenced using PCR primers as sequencing primers and the 
Applied Biosystems Prism BigDye terminator cycle sequencing 
reaction kit. The products were evaluated using an Applied 
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M 1 2 3 4 




655 bp 



Figure 5. Agarose gel electrophoresis of PCR product with the 
breakpoint primers F7 and R7. There is a 655-bp band in the 
patient lanes (1 and 2) but no product in the lanes of the unaffected 
family member (3 and 4). M: marker. 
doi:10.1371/journal.pone.0094201.g005 



Biosystems 3100 DNA sequencer. In addition, we genotyped the 
flanking polymorphic microsatellite markers of BMPRIB and 
GDF5 and performed a two-point linkage analysis using the 
LINKAGE 5.2 software program assuming a genetic model 
featuring autosomal dominance, a disease-allele frequency of 
0.0001, an evenly shared allele frequency, a zero phenocopy rate, 
no sex difference, and full penetrance [56] . 



Genome-wide scan and linkage analyses 

We performed a genome-wide scan of the family using a set of 
403 polymorphic microsatellite markers (ABI Prism Linkage 
Mapping Set version 2.5) with an average marker density of 10 
centimorgans according to the Marshfield genetic linkage map. 
Additional markers were used for the fine mapping. The products 
of the PCR assays with fluorescently labeled primers were 
analyzed by automated capillary genotyping on MegaBACE 
1000 (Amersham Pharmacia Biotech). A two-point linkage analysis 
was performed using the LINKAGE 5.2 software program 
[18,56]. The haplotypes were reconstructed using the CYRILLIC 
version 2.1 software. 

Mutation analyses of BMP2 

Following the procedures described above, we amplified and 
sequenced the BMP2 gene (i.e., all exons, introns, promoter 
region, 5'UTR, and 3'UTR). The evolutionary conserved regions 
that were sequenced were Chr20: 6,811,205-6,813,071. The 
following evolutionary conserved sequences were identified using 
the pubhcly available web-browser program ECR [57]. The PCR 
products were analyzed and sequenced by standard protocols as 
described above. 

Duplication analyses of BMP2 

A CGH array of human chromosome 20 was performed using a 
custom Roche array (NimbleGen Human CGH 385K Chromosome 
20 Tiling Array). The labeling and hybridization were performed 
following the protocols provided by die manufacturers (NimbleGen 
Hybridization System 4). The array was analyzed using the 
NimbleGen MS 200 Microarray Scanner and NimbleScan software. 
A graphical overview was obtained using the SignalMap software 
(vl. 9.0.05, NimbleGen Systems Inc). Quantitative real-time PCR (Q; 
PCR) with five primer pairs (P1-P5) was performed using a ViiATM 7 
Real-Time PCR system and the ViiATM 7 software 1.0 (Applied 
Biosystems) to confirm the duplication region. PCR with the primers 



Chr20:6,808,129 



Chr20:6,696,745-6,708,910 



5i>Id> 



Chr20:6,814,024 



Chr20:6,808,477 



5.5 kb 



BMP2 



Chr20:6,809,218 Chr20:6,813,888 
i 4,671bp 



Tel 



Cen 









4.5 kb 


^ 




• 


4.6 kb 







Chr20:6.809,382 



Chr20:6,814,044 



Figure 6. The chromosomal positions of the three previously reported duplications (Chr20: 6,808,129-6,814,024; Chr20: 6,808,477- 
6,814,024; and Chr20: 6,809,382-6,814,044) are indicated by the horizontal arrows. Comparison with the reported duplications shows 
that our duplication (Chr20: 6,809,218-6,813,888) overlaps with them but has different ends. Furthermore, the causal sequences of BDA2 could be 
narrowed down to 4,507 bp (Chr20: 6,809,382-6,813,888). 
doi:1 0.1 371/journal.pone.0094201 .g006 
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F7 and R7 was performed to amplify the breakpoint region, and the 
products were sequenced as described above. The abovementioned 
primers were designed by Su et al. [23] . The sequences were compared 
with the data from the UCSC Genome Browser (http:/ /genome.ucsc. 
edu). 
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